Constraint propagation in (N + l)-dimensional space-time 
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Higher dimensional space-time models provide us an alternative interpretation of nature, and 
give us different dynamical aspects than the traditional four- dimensional space-time models. 
Motivated by such recent interests, especially for future numerical research of higher-dimensional 
space-time, we study the dimensional dependence of constraint propagation behavior. The 
TV + 1 Arnowitt-Deser-Misner evolution equation has matter terms which depend on TV, but the 
constraints and constraint propagation equations remain the same. This indicates that there would 
be problems with accuracy and stability when we directly apply the TV + 1 ADM formulation to 
numerical simulations as we have experienced in four-dimensional cases. However, we also conclude 
that previous efforts in re-formulating the Einstein equations can be applied if they are based on 
constraint propagation analysis. 
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I. INTRODUCTION 

Higher dimensional space-time models have been in- 
vestigated from many viewpoints in physics. Current re- 
search interests come from brane-world models that try 
to solve the hierarchical problem in the unified theory 
(e.g. [1, 2]). Since these models can be probed by future 
Large Hadron Collider experiments, a lot of research is 
being undertaken. Even apart from such brane-world 
models, many new physical results in higher dimensional 
general relativity are reported. Although we do not have 
space to list them all, we mention the discoveries of the 
black-hole solutions in five-dimensional space-time (e.g. 
[3]) that violate the traditional black-hole no-hair conjec- 
ture, the possibility of new stable configurations of black- 
string models (e.g. [4, 5]), and the modified version of 
cosmic hoop conjecture (e.g. [6]). 

In order to investigate such topics, especially their dy- 
namical and nonlinear behavior, numerical simulations 
are necessary. Numerical relativity is promising research 
field, but it is also true that we have not yet obtain the 
recipe to perform long-term stable and accurate dynam- 
ical evolution. Many trial simulations of binary compact 
objects have revealed that the mathematically equivalent 
sets of evolution equations are showing different numer- 
ical stability in the free evolution schemes. Current re- 
search target in numerical relativity is to find out better 
reformulation of the Einstein equation (see reviews, e.g. 
[7, 8, 9]). 
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In this Report, we study the dimensional dependence of 
constraint propagation in the standard Arnowitt-Deser- 
Misner (ADM) formulation of the Einstein equation 
(space-time decomposition) [10, 11]. Reader might think 
that starting with the ADM equation is old-fashioned 
since recent large-scale numerical simulations are not us- 
ing the ADM equation due to its stability problem. How- 
ever, we still think that ADM is the starting formulation 
for analyzing the dynamical behavior both analytically 
and numerically. The plenty of re-formulations of the 
Einstein equations have been proposed in a last decade. 
Most of them are starting from the ADM variables. The 
practical advantages of such re-formulations are exten- 
sively under investigation by many groups now, but, in 
our viewpoint, the essential improvements of them can 
be explained in a unified way via constraint propaga- 
tion equations [8]. As we have shown in [12, 13], the 
stability problem of ADM can be controlled by adjust- 
ing constraints appropriately to evolution equations, and 
that key idea also works in other formulations [14, 15]. 
Therefore the analysis of the ADM equation is still es- 
sential. 

The idea of constraint propagation (originally reported 
in [16, 17]) is a useful tool for calibrating the Einstein 
equations for numerical simulations. The modifications 
to the evolution equations change the property of the as- 
sociated constraint propagation, and several particular 
adjustments to evolution equations using constraints are 
expected to diminish the constraint violating modes. We 
proposed to apply eigenvalue analysis to constraint prop- 
agation equations and to judge the property of the con- 
straint violation. The proposed adjusted equations have 
been confirmed as showing better stability than before 
by numerical experiments (e.g. [18, 19]). The purpose of 
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this report is to show this idea is also applicable to all 
higher dimensional cases. 

II. N + 1-DIMENSIONAL ADM EQUATIONS 

We start from the (N + l)-dimensional Einstein equa- 
tion, 

G„ v = 7^ - Ig^n + g„ v A = 8nT pu , (2.1) 

and decompose it into TV-dimensional space plus time, 
using the projection operator _L£, 

7£ = 5S+n"n„ = -L£. (2.2) 

where is a unit normal vector of the spacelike hyper- 
surface S, and we write the metric components, 

ds 2 = -a 2 dt 2 + 7 y {dx l + /3 i dt)(dx j + (3 ] dt), (2.3) 

where jij expresses iV-dimensional intrinsic metric, and 
a and f3 l the lapse and shift function, respectively. 
(Greek indices proceed 0,1, - • • ,N, while Latin indices 
proceed 1, • • • , N.) 

The projections of the Einstein equation are the fol- 
lowing three: 

G pu n»n v = toT^nPn" = 8np H , (2.4) 
G^n»±» p = torT^n* !?„ = -SnJ p , (2.5) 
G^^IZ = 8nT^^ p ^=8TTS prT , (2.6) 

where we defined 

Tp V = Phu^Uv + Jnn v + J v n^ + S^, (2.7) 

which gives T = —pa + S^g. 

To express the decomposition, we introduce the extrin- 
sic curvature as 

Kij = --L?-L?V„n„ 

= (-dtlij + Djfti + Dify) , (2.8) 

where V and Di is the covariant differentiation with re- 
spect to g^ v and 7^-, respectively. 

The projection of the Einstein equation onto £ is given 
using the Gauss equation, 

Wff jkl = - K\K 3l + K\K jk , 

(2.9) 



and the Codazzi equation, 

DjK^ - DiK = -K pa n a Lf , (2.10) 

where K — K l i. For later convenience, we contract (2.9) 
to 

(N) R = n + m^n" -K 2 + K ij K ij . (2.12) 
Eqs. (2.4) and (2.12) give the Hamiltonian constraint, 
Ch ~ 0, where 

C H = {G^-87TT^)nV 

= + k 2 - K^Kij) - 8np H - A,(2.13) 

while (2.5) and (2.10) give the momentum constraint, 
Cmi ~ 0, where 

C Ml = (G^-8^K^ 

= DjKj - DiK - 8irJ z . (2.14) 

Both (2.13) and (2.14) have the same expression as those 
of the four-dimensional version. 

The evolution equation for 7^ is obtained from (2.8), 
which is again the same expression as the four- dimen- 
sional version. 

The evolution equation of K i3 is obtained also from 
(2.5). The contraction of (2.1) gives 

Kii = ^ [Sij - 7^T7 y T ) - y-^7 y A{2.15) 

where we used g p „g^ u = N + 1. A straightforward calcu- 
lation of = d p Y% - 8 3 T% + Y^ - r^Tf p , where 

T^ p is the Christoffel symbol, gives 

all^pXnP = {d t K l3 ) + (DjDia) - (3 k (D k K l j) 
-(Dj(3 k )K ik - (Dif3 k )K kj + aK kl K k . (2.16) 

Substituting (2.15) and (2.16) into (2.11), we obtain 
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d t Kij = a {N) Rij + aKKij - 2aK t j K u - D t D a 
+p k (D k K l0 ) + (D (3 k )K lk + {B t fi k )K k , 
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-87TQ [ Sij - — 



2a 
N - 1 



(2.17) 



that is, only the matter and cosmological constant related 
terms depend on the dimension N . 

If we have matter, we need to evolve the matter thems 
together with metric. The evolution equations for mat- 
ter terms can be derived from the conservation equation, 
VT^j, = 0. In the next section, we will discuss them 
together with the constraint propagation equations. 



derivation via the Bianchi identity [17] may be the easi- 
est. 

In general, we write a N + 1-dimensional symmetric 
tensor which obeys the Bianchi identity, ^7 U S^ V = 0. 
Let us express by decomposing as 



III. N + 1-DIMENSIONAL CONSTRAINT 
PROPAGATION 

The constraint propagation equation, dt(CH ,Cmi) t , 
can be derived in many ways, and among them the 



S ilv =Xn jl n v + Y jl n v +Y v n il + Z jlv . (3.1) 



The normal and spatial projections of become 



n"V5^ = - V(W) - V*% + Y„n»VX ~ ^n v {V v n^) - A(V^) - n^X), (3.2) 

WVS^ = v^z v + r l (v^n M )+y M (v^)+x(v^K+n Al (v^r l ), (3.3) 

where we used V, while [17] uses different operator. For convenience, we rewrite them 

an"V5„„ = ~{d t X) + aKX + F(djX) - 07* (^V-) + a(^7mn)(7 m V j - (l/2) 7 m V j )^ 

-2 7 im (a m a)y, + aK^Zij, (3.4) 
-ahtVS^ = -(dtYi) - (d ia )X + aKYi + P(djYi) + 1 km {d i p m )Y k - & 7^ \d llpj )Y k 

-a^ k {8 k Z l0 ) - {dja)Zj + {l/2)a{d iljk )Z k i + ai mk {d mlkj )ZJ - (l/2)a 1 mk (d nmk )Z l \(i.b) 

respectively. ^ 



If we substitute (S M „ , X, Y { , Z i3 ■) = (T^, p H , Ji, <%) 
into (3.4) and (3.5) and assume V^T^^ = 0, then 
we obtain the matter evolution equations, dtpn and 
d t Ji- If we substitute (5 M „, X, Y i} Zij) = (G M „ - 

with K^const. and assume 



87rT M1/ , Ch, Cmu KjijCn) 



V^(G 



8irT, 



, then we obtain the constraint 
propagation equations, d t Cu and dtCui- [The parameter 
k corresponds to adding a term to (2.17), +(k — 1)Ch- ] 

This derivation does not depend on the dimension N 
at all. Therefore the evolution equations both for the 
matter and constraints remain the same with those in 
the traditional four dimensional version. 

The constraints include the extrinsic curvature terms, 
and the evolution equation of Kij changes due to N as we 
saw in (2.17). Interestingly, however, such changes will 
be cancelled out and the resultant constraint propagation 
equations remain the same. 

This means that a series of constraint propagation 
analyses can be directly applied to higher dimensional 
space-time. That is, the standard ADM evolution 
equations are likely to fail for long-term stable simula- 



tions. However, previously proposed adjustment tech- 
niques (e.g. [12, 13]) are also effective. 

For example, constraint amplification factors (i.e. 
the eigenvalues of constraint propagation matrix) in 
Schwarzschild space-time [eq. (47) in [13]] are 
(0,0, ±/(r)) for four-dimensional standard ADM evo- 
lution equations, where f(r) is a complex- valued func- 
tion. In the five-dimensional Schwarzschild or black- 
string case, they become simply (0, 0, 0, ±/(r)). 



IV. REMARKS 

Motivated by the recent interests in higher dimensional 
space-time, we checked the constraint propagation equa- 
tions based on the N + 1 ADM scheme. The evolu- 
tion equation has matter terms which depend on N, but 
we show the constraint propagations remain the same 
as those in the four-dimensional ones. This indicates 
that there would be problems with accuracy and stabil- 
ity when we directly apply the N + 1 ADM formulation 
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to numerical simulations as we have experienced in four- 
dimensional cases. However, we also conclude that pre- 
vious efforts in re-formulating the Einstein equations can 
be applied if they are based on constraint propagation 
analysis. The generality holds for other systems when 
their constraints are writtin in the form of (3.1). 

Since we only used the Bianchi identity in the core 
discussion, the assertion is also applicable to branc- 
world models. In the context of the Randall-Sundrum 
brane- world models [2], people study the modified four- 
dimensional Einstein equations [20], which are derived 
from five-dimensional Einstein equations with a thin- 
shell (3-brane) approximation. The terms there addi- 
tional to the standard ADM (see eq.(17) in [20]) in- 
clude extrinsic curvature (due to shell- normal vector), 
cosmological constant (s), and five-dimensional Weyl cur- 



vature. These terms, however, can be interpreted as a 
single stress-energy tensor which obeys the Bianchi iden- 
tity. Therefore the properties of the constraint propa- 
gation equations are the same as the above (from the 
five-dimensional space-time viewpoint). Our proposals 
for the adjustments [12, 13] are also valid in brane-world 
models. 

We hope this short report helps numerical relativists 
for developing their future simulations. 

Acknowledgements HS is supported by the special 
postdoctoral researchers' program at RIKEN. This work 
was supported partially by the Grant-in- Aid for Scientific 
Research Fund of the Japan Society of the Promotion of 
Science, No. 14740179 and by Waseda University Grant 
for Special Research Projects Individual Research 2003A- 
870. 



[1] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. 
Lett. B 429, 263, (1998); ibid. 436, 257 (1998). 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 
(1999); ibid. 83, 4690 (1999). 

[3] G. T. Horowitz and A. Strominger, Nucl. Phys. B. 360 
197 (1991); R. Emparan and H. S. Reall, Phys. Rev. Lett. 
88 101101 (2002). 

[4] R. Gregory and R. Laflamme, Phys. Rev. Lett. 70, 2837 
(1993); G. T. Horowitz and K. Maeda, Phys. Rev. Lett. 
87, 131301 (2001). 

[5] T. Wiseman, Class. Quant. Grav. 20, 1177 (2003); M. 
Choptuik, L. Lehner, I. Olabarieta, R. Petryk, F. Preto- 
rius, and H. Villegas, Phys. Rev. D 68, 044001 (2003). 

[6] D. Ida and K. Nakao, Phys. Rev. D 66, 064026 (2002). 

[7] L. Lehner, Class. Quantum Grav. 18, R25 (2001). 

[8] H. Shinkai and G. Yoneda, in Progress in Astronomy and 
Astrophysics (Nova Science Publ., New York) to be pub- 
lished. The manuscript is available as gr-qc/0209111. 

[9] T. W. Baumgarte and S. L. Shapiro, Phys. Rep. 376, 41 
(2003). 

[10] R. Arnowitt, S. Deser and C.W. Misner, in Gravitation: 
An Introduction to Current Research, ed. by L.Witten, 
(Wiley, New York, 1962). 



[11] J.W. York, Jr., in Sources of Gravitational Radiation, 
ed. by L. Smarr, (Cambridge, 1979) ; L. Smarr and J.W. 
York, Jr., Phys. Rev. D17, 2529 (1978). 

[12] G. Yoneda and H. Shinkai, Phys. Rev. D 63, 124019 

(2001) . 

[13] H. Shinkai and G. Yoneda, Class. Quantum Grav. 19, 
1027 (2002). 

[14] H. Shinkai and G. Yoneda, Class. Quant. Grav. 17, 4799 

(2000) . G. Yoneda and H. Shinkai, Class. Quant. Grav. 
18, 441 (2001). 

[15] G. Yoneda and H. Shinkai, Phys. Rev. D 66, 124003 

(2002) ; Class. Quant. Grav. 20, L31 (2003). 
[16] S. Detweiler, Phys. Rev. D 35, 1095 (1987). 
[17] S. Frittelli, Phys. Rev. D 55, 5992 (1997). 

[18] B. Kelly, P. Laguna, K. Lockitch, J. Pullin, E. Schnetter, 
D. Shoemaker and M. Tiglio, Phys. Rev. D 64, 084013 

(2001) . 

[19] H-J. Yo, T.W. Baumgarte and S.L. Shapiro, Phys. Rev. 

D 66, 084026 (2002). 
[20] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D 

62, 024012 (2000). 



